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^ ; Abstract 

C ' The preparation of mesoscopic states of the radiation and matter fields through atom-field in- 

teractions has been achieved in recent years and employed for a range of striking applications in 
quantum optics. Here we present a technique for the preparation and control of a cavity mode 
which, besides interacting with a two-level atom, is simultaneously submitted to linear and para- 
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■ metric amplification processes. The role of the amplification-controlling fields in the achievement 

: 

O ■ of i^sal mesoscopic states, is to produce highly-squeezed field states and, consequently, to increase 

m 

I both: i) the distance in phase space between the components of the prepared superpositions and 

pi • ii) the mean photon number of such superpositions. When submitting the squeezed superposition 

. states to the action of similarly squeezed reservoirs, we demonstrate that under specific conditions 

^ , the decoherence time of the states becomes independent of both the distance in phase space be- 



tween their components and their mean photon number. An explanation is presented to support 



^ . this remarkable result, together with a discussion on the experimental implementation of our pro- 

H : 

. . , posal. We also show how to produce number states with fidelities higher than those derived as 

circular states. 
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I. INTRODUCTION 



The successful manipulation of atom-field interactions in cavity quantum electrodynamics 
(QED) and trapped ions is a great achievement of present-day physics which has encouraged 
outstanding theoretical proposals and experimental implementations. As high-Q cavities [l| 
and ionic traps P| have permitted the preparation of coherent-state superpositions of the 
form = (jae*'^) + lae"**^)) /"\/2, with mean numbers of photon and phonon quanta 
\a\ ~ 10, mesoscopic quantum coherence has been investigated. In the cavity QED do- 
main, the progressive decoherence of mesoscopic superpositions involving radiation fields 



1| and the re- 
3] . Moreover, 



with classically distinct phases was observed through atom-field interaction 
versible decoherence of such a mesoscopic-field state has been conjectured 
the generation and detection of Fock-states of the radiation field was demonstrated experi- 
mentally Q] and the Rabi oscillation of circular Rydberg atoms in the vacuum and in small 
coherent fields in a high-Q cavity was measured ^|, revealing the quantum nature of the 
radiation field 0|. 

Parallel to the achievements in cavity QED, the mastery of techniques to manipulate 
electronic and motional states of trapped ions with classical fields has enabled the control of 
fundamental quantum phenomena at a level that seems to herald a new phase in technology. 
The operation of a two-bit controlled-NOT quantum logic gate was demonstrated by storing 
the two quantum bits in the internal and external degrees of freedom of a single trapped 
ion [3]. A "Schrodinger cat" superposition of spatially-separated coherent harmonic states 
was generated as well as other nonclassical states at the single atom level j^. The 
reconstruction of the density matrices and the Wigner functions of various quantum states 
of motion of a harmonically bound ion was also reported j^. 

Besides the atom-field interaction in cavity QED and trapped ions, the preparation of ref- 
erence travelling-field states needed to measure the properties of signal travelling fields 10| 
has been suggested based on optical linear and nonlinear devices The techniques 
developed over the last decades for the process of parametric up- and down-conversion have 
enabled great advances in the domain of travelling waves. The production of true entangle- 
ment by type-II noncoUinear phase-matching in parametric down-conversion was employed 
to demonstrate a violation of Bell's inequality with two-photon fringe visibilities in excess of 
97% 13]. Three-photon Greenberger-Horne-Zeilinger entanglement has also been observed 
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Ql and it is worth stressing that the experimental implementations of teleportation have 
been achieved with travelling wave techniques as these provide the facilities for prepa- 
ration of the state to be teleported, for the quantum channel and for the accomplishment of 
the required Bell-type measurements. High-fidelity teleportation of superpositions [16] and 
entanglements of running-wave field states ^ have also been presented. 

As the techniques for generating nonclassical superposition states have been improved, 
the attention turns to a major problem that must be overcome in the contexts of quantum 
communication [3] and computation [19]: the decoherence of qua ntum states due to the 
inevitable coupling of the quantum systems to their environment [21] 
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and even due to 



fluctuation in the interaction parameters required to prepare a coherent superposition 
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l24l ]. Schemes for inhibiting decoherence by engineerin g th e reservoir have been developed for 
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Measurements of the decoherence 



trapped ions j25|, |26[ and atomic two-level systems 
of superposed motional states of a trapped ion coupled to an engineered reservoir, where the 
coupling and the state of the environment are controllable [2^, have also been reported. The 
possibility of controlling the decoherence mechanism is crucial to the preparation of the long- 
lived macroscopic superposition states and entanglements of macroscopic objects required for 
the implementation of the potential applications of quantum mechanics in communication 
and computation [3^ . Beyond the efforts being made to engineer mesoscopic superpositions 
and entanglements with photon and phonon quanta, referred to above, the possibility of 
engineering such mesoscopic states with massive objects has been and is being pursued. 
Recently, correlations (on average) of collective spin states of two macroscopic objects, each 
consisting of a caesium gas sample with about 10^^ atoms, was demonstrated experimentally. 
In Ref. j3l| an experimental scenario designed to reduce dramatically the decoherence rate of 
a quantum superposition of Bose-Einstein condensates is outlined. This is also our concern 
in the present work, focusing on the preparation of long-lived states of the radiation field in 
cavity QED. 

Methods for manipulating cavity-field states through atoms driven by external fields 
constitute an important means of attaining arbitrary control of a quantum field. Although 
the time evolution of a field state under linear and parametric amplifications has been a major 
concern in quantum optics for generating squeezed states and investigating their properties 



33, 



of classical driving fields have barely been considered for quantum states engineering 



purposes. Here we present a proposal for achieving long-lived mesoscopic superposition 
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states of the radiation field in cavity QED wliich relies on two basic requirements: parametric 
amplification and an engineered squeezed-vacuum reservoir for cavity-field states (we note 
that the required engineered reservoir - resulting from the standard vacuum for cavity modes 
plus additional interactions - must be an optimum squeezed- vacuum reservoir). In addition, 
our technique can be emplo yed to prepare number states with fidelities higher than those 
generated as circular states |3q . 

Our proposal considers the dispersive interaction of a two-level atom with a cavity field 
which is simultaneously under amplification processes. The parametric amplification is em- 
ployed to achieve a high degree of squeezing and excitation of what we actually want to be 
a mesoscopic superposition state. We show that the prepared squeezed-mesoscopic state, 
under the action of a similarly squeezed reservoir, exhibits a decoherence time orders of 
magnitude longer than those of non-squeezed cavity-field states subjected to the influence 
of i) a squeezed reservoir and ii) a non-squeezed reservoir. In fact, the computed decoher- 
ence time turns out to be independent of both the average photon number and the distance 
in phase space between the centers of the quasi-probability distribution of the individual 
states composing the prepared superposition. The decoherence time depends only on the 
excitation of the initial coherent state injected into the cavity previous to the preparation 
of the squeezed superposition. This remarkable result follows when the direction of squeez- 
ing of the superposition state is perpendicular to that of the reservoir modes. Under this 
condition, the entanglement between the prepared state and the modes of the reservoir is 
minimized and so the noise injected from the reservoir into the prepared cavity mode is 
minimal, making it a long-lived superposition state. 

We finally stress that a scheme has been presented in Ref. js^ for the implementation 
of the parametric amplification of an arbitrary radiation-field state previously prepared in 
a high-Q cavity. As squeezed l igh t is mainly supplied by nonlinear optical media as run- 
ning waves (through backward [22| or forward four-wave mixing and parametric down- 
conversion [221), standing squeezed fields in high-Q cavities or ion traps can be generated 
through atom-field interaction [40]. Although considerable space has been devoted in the 
literature to the squeezing process in the Jaynes-Cumming model, the issue of squeezing any 
desired prepared cavity-field state |\E'), i.e., the accomplishment of the operation S{() |\E') 
in cavity QED {( standing for a set of group parameters) has not been addressed. Engi- 
neering such an operation was the subject of Ref. 0]; it is achieved through the dispersive 
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interactions of a three-level atom simultaneously with a classical driving field and a cavity 
mode whose prepared state we wish to squeeze. In short, the dispersive interaction of the 
cavity mode with a driven atom produces the desired operation S{() |\E'). Since linear am- 



plification is easily accomplished in cavity QED |41|, |42| , the scheme in Ref . 36| contributes 
crucially for the experimental feasibility of the present proposal for preparation and control 
of long-lived cavity-field state through atom-driving field interaction. 



II. ATOM-DRIVEN FIELD INTERACTION 



The proposed configuration for engineering driven- cavity-field states, depicted in Fig. 1, 
consists of a two-level Rydberg atom A which crosses a Ramsey-type arrangement, i.e., a 
high-Q micromaser cavity C located between two Ramsey zones Ri and R2. After interacting 
with this arrangement, the atom is counted by detection chambers D2 and Di (for ionizing 
the excited |1) and ground |0) states, respectively), projecting the cavity-field in the desired 
state. The transition of the two-level atom A from excited to ground state is far from 
resonant with the cavity mode frequency, allowing a dispersive atom-field interaction to 
occur. In addition to the dispersive interaction with the two-level atom, the cavity mode is 
simultaneously submitted to linear and parametric amplifications (both represented in Fig. 
1 by the source S) so that the Hamiltonian of our model (for ^ = 1) is given by 

H = u^a^a + ^a. + xa^aa. + Ti^.,..^..^. (1) 

where = |1)(1| - |0)(0|, and a are, respectively, the creation and annihilation operators 
for the cavity mode of frequency uj which lies between the two atomic energy levels, which are 
separated by Uq, such that the detuning 6 = \u! — ujq\ is large enough to enable only virtual 
transitions to occur between the states |0) and The atom-field coupling parameter inside 
the cavity is x = ^^/^ j where Q is the Rabi frequency. The expression for the atom-field 
dispersive interaction on the right-hand side (rhs) of Eq. is valid under the assumption 
that Q'^n <^ 5^ + 7^, where n is a characteristic photon number and 7 is the spontaneous- 
emission rate j^]. We suppose, for simplicity, that the atom-field coupling is turned on (off) 
suddenly at the instant the atom enters (leaves) the cavity region, such that x = when 
the atom is outside the cavity. 

We consider the atom, prepared at time to by the Ramsey zone Ri in a |0),|1) super- 
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position, to reach C at time ti and leave it at t2- The hnear and parametric pumping are 
assumed to be turned on also at to and turned off at a convenient time t > t2- Finally, 
the action of the classical amplification mechanism on the cavity mode is described by the 
Hamiltonian 

nampUficaUon = C{t)a^' + C{t)a + + C{t)a, (2) 

where the time-dependent (TD) functions ({t) and ^(t) allow the parametric and linear am- 
plifications, respectively. It is well understood that for specific values of these TD functions 
the eigenstates of the Schrodinger equation may squeeze the variance in one of the cavity 
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modes' two quadrature phases [221 1 

The Schrodinger state vector associated with Hamiltonian can be written using 

where |$£ (t)) = / ^-^e (ct^'t) \ce), i = 0,1, the complex quantity a standing for the eigen- 
values of a, and Ae{a,t) = (a, £ |\E' (t) ) are the expansion coefficients for \^e(t)) in the 
coherent-state basis, {!«)}. Using the orthogonality of the atomic states and Eqs. dH) and 
dni) we obtain the uncoupled TD Schrodinger equations: 

z||<l>,(t)) = H,|$,(t)), (4) 

He = uJeit)a^a + C(t)a"f' + Cit)a + ^{t)a^ + C{t)a, (5) 

with ujeit) = uj — (— l)'^x0(t — ti)0(t2 — t) ■ Note that the problem has been reduced to 
that of a cavity field, under parametric and linear pumping, whose frequency uj is shifted by 
~X when interacting with the atomic state (1), during the time interval t = t2 — ti. 

Solving Eq. (0)) we obtain, from an initial state of the cavity mode at time t,, |$^()f:j)), 
the evolved state 

\^,{t)) = V,{t,U)\^,{U)) , (6) 

which defines the evolution operator we are looking for. Evidently, the evolution operators 
U£(ti, to) and U^(t, ^2)5 giving the evolution of the state vector of the radiation field while the 
atom is outside the cavity, do not depend on the state of the two-level atom, the label ^ being 
unnecessary. However, the operator U^(t2,ti), which gives the evolution of the cavity-field 
state during its interaction with the atom, does depend on the atomic state and differs from 
the operators U(ti,to) and 0(^,^2) only by the shifted frequency uje,{t). 



III. SOLVING THE SCHRODINGER EQUATION VIA TIME-DEPENDENT IN- 
VARIANTS 



The Hamiltonian in Eq. (0) has been inve stig ated in the search for squeezed states 



47[ and TD invariants 45[ have been 



of the radiation field. Group-theory methods 
used in attempts to solve this TD quadratic Hamiltonian, which may represent a charged 
particle subjected to a harmonic motion, immersed in a TD uniform magnetic field, a single 
mode photon field travelling through a squeezing medium or, as in the present situation, 
a cavity mode with shifted frequency under linear and parametric amplification. In the 
present work, we make use of the TD invariants of Lewis and Riesenfeld 4^ to solve the 
Schrodinger equation (0]), following the reasoning in Ref. 45|: instead of proposing an 
invariant associated with the Hamiltonian (jS)), we first perform a unitary transformation on 
Eq. ^ in order to reduce it to a form which already has a known associated invariant. 
Thus, under a unitary transformation represented by the operator S{e£) [e^ standing for a 
set of TD group parameters which may also depend on the atomic state i) , we obtain from 
Eq. dH) 



.d 

''di 



|$f(t)> = Hf |$f (t)> 



where the transformed Hamiltonian and wave vector are given by 



dt 



|$f (t)> = St(5,)|$,(t)). 



(7) 

(8a) 
(8b) 



In what follows we employ two theorems to obtain the solution of the TD Schrodinger 
equation a) a theorem expounded in j3l asserts that if Ie{t) is an invariant associ- 
ated with (i.e., dU{t)/dt = dU/dt + i\H^,h{t)\ = 0), then the transformed operator 
//(t) = S''{6£)Ie(t)S{6£) will be an invariant associated with ?if ; b) on the other hand, from 
Lewis and Riesenfeld's well-known theorem 4^ , it follows that a solution of the Schrodinger 
equation is an eigenstate of the Hermitian invariant /^(t) multiplied by a TD phase factor. 
It follows from a) and b) that the solutions of Eq. Q are given by 



\^i,m{t)) = S{ee)\<^l^{t)) = S{ee)e''^l^^'^ \m,t)s , m = 0,l,2,.... 



(9) 
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where |m, t)^ is the eigenstate of the invariant j49[ and the Lewis and Riesenfeld phase |48l | 
obeys 

^Imit) = 1^ dt's (m, t'\ (^t^ - \m, t'), . (10) 

It is straightforward to verify that under the unitary transformation carried out by the 
operator S{e£) the TD phase is invariant: 0fm(t) = 4>£^m(t). 

A. The transformed Hamiltonian 



Next, we associate the unitary transformation with the squeeze operator S{ee) 



exp 



I yega'' — e^a'^j , where the complex TD function eg = rf {t) e**^^*^*) includes the squeeze 
parameters ri{t) and (pe{t). ri{t) is associated with a squeeze factor, while (pi{t) defines the 
squeezing direction in phase space. Moreover, the TD parameters for the parametric and 
linear amplification processes are written as C,{t) = ^(t) e*''*^*^ and ^(t) = x(t) e*^*^*\ respec- 
tively. The squeeze parameters (^^(t), ipe{t)), the amplification amplitudes (K(t),x(t)) and 
frequencies {r]{t),tu{t)) are all real TD functions. With the above assumptions and after a 
lengthy calculation, the transformed Hamiltonian becomes 

= ni{t)a^a + Ae{t)a^ + k}{t)a + Fiif), (11) 

provided that its TD coefficients satisfy 

Vt^it) = oJiit) + 2K{t) tanh r^(t) cos (r/(t) - ifiit)) , (12a) 
A^(t) = i{t) coshr^(t) + C{t) e''^'^*^ sinhr^(t), (12b) 
Flit) = K{t) tanhr£(t) cos {7]{t) - ^e{t)) , (12c) 

while the squeeze parameters r^it) and feit) are determined by solving the coupled differ- 
ential equations 

reit) = 2K{t) sin {r]{t) - i^,{t)) , (13a) 
(p^{t) = -2uJi{t) - AK{t) coth (2r£(t)) cos (r/(t) - (pe{t)) . (13b) 

It is evident from these relations that the TD group parameters se{t), defining the unitary 
operator S{ee), depend on the atomic state £, as assumed from the beginning. We finally 
mention that we have associated the unitary transformation with the squeeze operator since 
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the parametric amphfication described by Hamiltonian (j2I) actually squeeze the cavity-field 
state. In fact, the TD parameter ({t) allowing the parametric amplification in Eq. Q is 
connected to the squeeze parameters {ri(t), fiit)) as expressed by Eqs. (|13ap and ()13b|) . 



B. The evolution operators 

With the Hamiltonian (jlip at hand we return to the solution of the Schrodinger equation 
(|7j). The application of the invariant method leads to the wave vector 4^ 

\^Lit)) = e"t'^'-^''>Dm)]\m) , m = 0,l,2,..., (14) 

where \m) is the number state and D [9i{t)] = exp \9i(t)a'^ — ^|(t)a] is the displacement 
operator, 6i{t) being a solution to the equation i9i(t) = Qi(t)9e{t) + Ai{t), given by 

e^it) = e-^^^W (e,{U) - I f A,{t') e^'^^(*') dt') , (15) 



with Peii) = /j* ^t{t')dt' . We note that 6'^(to) describes the initial cavity-field state which 
will be assumed to be a coherent state \a), the subscript C. being purely formal. From the 
substitution of Hamiltonian (fTTj) into the Lewis and Riesenfeld phase, defined in Eq. (fTUI). 
we obtain 

= - £ \^niQeit') + ^ [Al{t')ei{t') + Ae{t')em] + ri{t')'^dt'. (16) 

Therefore, the solutions of the Schrodinger equation (0)), which form a complete set, can 
be written 

\^e,mit)) = S [seit)] \<^l^it)) = U,{t) \m) , (17) 

where 

U,{t) = T,{t)S [e,{t)] D [e,{t)] R [n,{t)] (18) 

is a unitary operator containing, in addition to the squeezed and the displacement operators, 
a global phase factor 

T,(t) = exp |-2 [Re [A\{t')e^{t')] + r^(t')] cit'j • (19) 

The rotation operator in phase space, derived from the TD Lewis and Riesenfeld phase 
factor, is given by 

R[ni{t)] = exp [-ia^a(3e{t)] , (20) 



Hence, for the solution of Schrodinger equation HJ we find 

oo oo 
I'^iit)) = Y,Crn = U,{t) Cm |m) = U ,{t)U\{t,) \^ ,{t,)) , (21) 

m=0 m=0 

which finally defines the evolution operators 

Veit,U) = Ue{t)Ul{U). (22) 

We note that for the initial time R[ni{0)] = R[0] = 1, D[9i{0)] = D [a], S[ee{0)] = 
S [0] = 1, and T^(0) = 1. 

IV. EVOLUTION OF THE ATOM-FIELD STATE 

Let us assume that the micromaser cavity is prepared at time to in a single-mode coherent 
state la) by a monochromatic source, such that with m = in Eq. (fTTj) we have 9i{to) = a. 
Classical microwave fields are injected into the cavity and the amplitudes of these fields can 
be adjusted by varying the injection time. As mentioned above, the linear and parametric 
pumping are supposed to be turned on, also at to; the same time the atom is prepared by 
the Ramsey zone Ri in the superposition state cq |0) + ci |1). The combined atom-field state 
at the initial time to is, from Eq. ^ 

1^ (to)) = [e'"°*°/' Co |0) + e-'"«*«/' ci |1)] |a) . (23) 

In fact, with Ae{P,to) = (/?, £| (cq |0) + Ci |1)) |a;) it follows immediately that \^e{to)) = 
J^A,iP,to)\P) = c,\a). 

The evolution of the initial state (to)) to the time the atom reaches the cavity reads 

|^(ti))=U(ti,to)|^(to)). (24) 

Evidently, the evolution operators I[J(ti,to) and U(t,t2), which govern the dynamics of the 
cavity-field state while the atom is outside the cavity, do not depend on the state of the two- 
level atom. On the other hand, during the time interval r = t2~^i the atom spends inside the 
cavity the evolution of the entire system is dictated by the operator U^(t2, ti) = Ueit2)Ujiti). 
This depends on the atomic state i and differs from the operators U(ti,to) and U(t, ^2) by 
the shifted frequency u!i{t). Therefore, at the time the atom leaves the cavity, the state of 
the atom-field system is given by 

\^{t2)) = [e*"°*^/2co|0)Uo(t2,ti)+e-*"«*^/2ci|l)Ui(t2,ti)]U(ti,to)|«). (25) 
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After crossing the cavity, the atom evolves freely from t2 until the time it reaches the second 
Ramsey zone, i?2- During this time interval, the cavity mode continues to be pumped and 
the complete state of the system, evolving under the operator 0(^,^2), reads 

I* (t)) = U(t, h) [e^"«*/' Co |0) Uo(t2, ti) + e-''^«*/' c^ \ 1) lJ^{h, t,)] U(ti, to) \a) . (26) 

Next, the atom crosses the Ramsey zone R2, where a n/2 pulse is applied, leading the 
atom-field system to the entangled state 

|M/ (t)) = -L {[-e-«*/2coUo(t,to) +e-^-°*/2ciUi(t,to)] |0) 

where we have defined the operators 

U,(t,to) =U(t,t2)U,(t2,ti)U(ti,to). (28) 

Finally, measurement of the atomic state projects the cavity field into the "Schrodinger 
cat" -like state 

|vl/(t)) =A4 [±e^-°*/2^oUo(t,to)+e-*"«*/2ciUi(t,to)] |«) , (29) 

where the sign + or — occurs if the atom is detected in state |1) or |0), respectively, and 
A/± refers to the normalization factors. From Eq. it follows that, after measuring the 
atomic level used to generate the superposition state of the radiation field, it is possible 
to control this superposition by adjusting the TD amplification parameters K{t), x(t), ri{t), 
and zu{t). 

It is worth noting that expression can be manipulated, employing Eqs. (PHj) . (|22I), 
and Eq. |TH|). to give the simple form 

|vl/(t)) = A4 [±e^-«*/2coTo(t)5[eo(t)] \9o{t)) + e-'^^'^' c,T,{t)S [e^it)] m))] 

1 

= Af±J2ce{t)S[ee{t)]D[9e{t)] |0) 

1 

= Af±Y.c,{t)S[ee{t)]m)), (30) 

where ce{t) = ± (±)^e'^~^ iujot/2 QX^(t) and the amplitude of the coherent state \6£{t)) follows 
from Eq. |T3|). 
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A. Passing N atoms through the cavity 



Let us proceed to the construction of a cavity-field state by passing two or more atoms 
tfirough cavity C. It is easy to conclude from Eq. (j29|) that, after the passage of atoms 
through cavity C, each atom prepared in the state Co,a; |0) + Ci^k |1) by /c = 1, A^, we 
obtain the cavity-field state 

N 

l^^N (t)) = A4 n [± e""*^' co,fcUo,fc(t/,fc, + e-''-«*/2 ci,,Ui,,(t^,fc, t,,,)] \a) , (31) 

k=l 

where ti^^ stands for the time when the fc-th atom is prepared by Ri and stands for the 
time when the k-th atom is detected, assumed to be the same as ti^k+i- Therefore, we obtain 
^e,k{'tf,k,ti,k) = U(t/,fc = ti^k+i,t2,k)l^e{t2,k,ti,k)l^{ti,k,ti,k = tf,k~i)- After some manipulation 
and using Eq. (jHUj) . the state (jHTj) can be simplified to the form 

2 N 

^l,.../]V=l fc=l 

= Af±Y.^,it)S[Ekit)]\Mt)), (32) 

k=l 

where we have replaced by XlLi' i-^-, ^^(t) = r£,,...,f^(t) exp (iv?£i,.../^r(t)) and 

V. ANALYTICAL SOLUTIONS OF THE CHARACTERISTIC EQUATIONS 
(13A,13B) 

In this section we present some specific solutions of the characteristic equations (jl3all3b|l . 
following a more denned .e..e:. . Q. We ,„ve.,a.e .Ke ..ua..: w.e. 

mode I a) is resonant with the driving classical fields during the time the atom is out of the 
cavity: from to to ti and from t2 to t. The parametric amplifier is assumed to operate in a 
degenerate mode in which the signal and the idler frequencies coincide, producing a single- 
mode driving field. In the resonant regime this single-mode field has the same frequency lj 
as the cavity mode so that rjit) = —2ujt 33]. For the resonant linear amplifier it follows 
that ^{t) = —ujt. However, during the time that the atom is inside the cavity, from ti to t2, 
it pulls the mode frequency out of resonance with the classical driving fields, establishing a 
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dispersive regime of the amplification process. Thus, in what follows we derive the solutions 
of the coupled differential equations ()13all3bj) for the resonant and dispersive regimes. 



A. Resonant amplification 

We start with the solution of the characteristic equations ()13all3bp for the resonant 
amplification that operates while the atom is out of the cavity, from to to ti and ^2 to t, when 
so that uJe{t) = u, w{t) = —ujt and ri{t) = —2ujt. Defining uj = f{t), ip{t) = —2f{t) + g{t), 
and 7]{t) = -2/(t) + h{t), Eqs. (I13all3b|l become 

r{t) = -2n{t) sin {g{t) - h{t)) , (33a) 
g(t) = -4K(t) coth (2r(t)) cos {g{t) - h{t)) . (33b) 

Assuming that h{t) = h is a constant, the time dependence is eliminated from Eqs. ()33ap . 
()33b|) and we are left with the first-order differential equation 

dv 1 

— = - tanh(2r) tan(^ - h). (34) 
ag Z 

After integrating Eq. (|34j) we obtain the constant of motion 

cos {(fit) - ri{t)) sinh (2r(t)) = (35) 

with Ci depending on the initial values r(tj), ^{U), and ^^(tj), where z = 0,2. Thus, the 
solutions of Eqs. ()33a|l and ()33b|l . which apply under the condition cosh^ (2?"(t)) > 1 + , 
are given by 



cosh (2r(t)) = V 1 + Cf cosh 



cosh-(^^^^^^|±.(M. 



(36a) 



C 

cos {^{t) - r]{t)) = ^ ' (36b) 



'cosh^2r(t) - 1 
where 

u{t,ti) = 4 / K{t)dt. (37) 



Note that for = the cosh ^ term in Eq. ()36ap is null, and the signals ± become irrelevant. 
However, for ti = t2 we choose the sign that gives r(t) > 0. 
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B. Dispersive amplification 



Dispersive amplification occurs during the time the atom is inside the cavity, shifting the 
mode frequency uj hj x = so that uj£ = u!±x- Evidently, the amplification frequencies 

are unaffected by the passage of the atom, so that r]{t) = —2ujt and w{t) = —uot. Assuming 
that parameter k is time-independent and defining (f£{t)—r]{t) = fi{t) and 77(t)+ 2uji{t) = g^, 
Eqs. p3all3b|) become 

n{t) = 2K{t) sin {fe{t)), (38a) 
fiit) = -g, - Mt) coth (2r,(t)) cos im) . (38b) 

Since g^ = — (— 1)^2% is a constant, Eqs. (|38al38b|) can be solved by quadrature, leading to 
a constant of motion 

cosh (2r,(t)) + cos Mt) - r,{t)) sinh (2r,(t)) = Ci, (39) 

which now depends on the initial values r(ti), ^p{ti), and rjiti). Despite the assumption that 
the atom- field coupling is turned on (off) suddenly, these initial values must be computed 
from the solutions for the resonant amplification regime at time ti. With this procedure we 
obtain the solutions for the resonant amplification (r(ti), ^p{ti)) as the limit of those for the 
dispersive amplification {r(ti), f(ti)) as x — ^ 0. The parameter = — (— 1)^2k/x, defined 
for a constant amplification amplitude k, is an effective macroscopic coupling. Therefore, for 
the dispersive regime, we find three different solutions, depending on whether the coupling 
is strong (l^^l > 1), weak (l^^l < 1), or critical coupling (|^^| = 1). 
a) With strong coupling (|^£| > 1), we have the relations 

-g/i{t) 



cosh {2re{t)) 



(40a) 



/ /X Ci — cosh (2r At)) X 



where 



Mi) = T^^^l^u(t,ti) + ln 



2 m (^(^^^iHcf^ + Ci 1^, 



1^. 

the sign being chosen so that r{t) > 0. The function u{t,ti) is defined by Eq. (|?fj). 



(41) 
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b) For the weak coupling regime < 1), the TD squeeze parameters when Ci > 

— are given by Eq. ()40b|) and 



cosh(2r£(t)) 



We 



sm 



1^^ 



Cx 



u{t, ti) + arcsin 



(42) 



c) Finally, for critical coupling (|^^| = 1), the TD squeeze parameters are given by Eq. 
fl40bj) and the solution for r^(t) follows from the equation 



cosh(2r^(t)) 



1 

2C'i 



l + Cf+ (^VCi(2cosh(2r(ti)) -Ci) - 1 T Ciu{t, ti) 



(43) 



the parameter i being redundant. Note that in Eqs. (j4(Jaj) and (j42j) the parameter i is also 
unnecessary since the rhs of Eq. (j4Uaj) is an even function of Therefore, the squeezing 
factor r{t) does not depend on the atomic state, in contrast to the squeezing direction in 
phase space defined by ipe(t). 

/^From the above solutions for the resonant and dispersive amplifications it is straightfor- 
ward to obtain the behavior of the TD squeeze parameters from time to, when the classical 
driving fields are turned on simultaneously to the preparation of the atomic state by Ri, to 
any instant t after the atom-field interaction. The time t may be chosen to be before, after 
or in the course of the ionization detection of the atomic state: 

i) From to to ti, the squeeze parameters follow from Eqs. ()36aj) and ()36b|) . As mentioned 
above, such equations apply under the condition cosh^ (2?"(i)) > I + Cf, which is always 
satisfied for Cj = 0, a value following from the initial conditions r(to), f(to), and ri(to). In 
fact, for an initial coherent state injected into the cavity: r{to) = 0. Assuming the parameter 
K to be time- independent, together with Co = 0, Eqs. (j36aj) and ()36b|) lead to the simplified 
solutions 



r{t) = 2Kt, 

ip{t) = -2ujt + 7r/2. 



(44a) 
(44b) 



ii) From ti to ^2 we have three possible solutions for the squeeze parameters, depending 
on the coupling strength \%^e\. These solutions follow from the above- described expressions, 
in Eqs. fl40a|) . ()40b|) . (jl^ . and (jlHj) . given that the constant of motion Ci = cosh(4Kti), 
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computed from Eqs. (jSHI), ()44a|) . and ()44b|) with t = ti. It is straightforward to observe 
in these equations the well-known threshold in the behavior of the TD squeeze factor r{t), 
following from the quadratic TD Hamiltonian (0) [4^: r{t) increases monotonically for 
l^^l ^ 1; while for < 1 it oscillates periodically. In the present paper we are interested 
in the weak coupling regime, where the squeeze parameters follow from Eqs. ()40b|) and (|42j) . 
We note that for realistic physical parameters we achieve higher squeezing factor even in 
this regime. 

iii) From ^2 to t the squeeze parameters are again derived from Eqs. ()36a|) and ()36b|) . 
The constant of motion is computed from the initial conditions ^(^2), f{t2), and ri{t2), which 
depend on the strong, critical or weak coupling regimes. For weak coupling, I^P^I < 1, in 
which we are interested, the constant of motion in Eq. (j35|) . derived from Eqs. ()40b|l and 
(jl^ and depending on the atomic state, reads 

C„ ^ (45) 

VI. WIGNER FUNCTIONS AND FLUCTUATIONS OF THE QUADRATURES 

Now we analyze the states (jHUj) and ((221) , projected into the cavity after the detec- 
tion of one or several atoms, respectively, and especially control of these states through 
the amplification parameters. From here on we assume that the atom is detected in 
excited state so that = ^f+J2]=oCi{t)S[£i{t)]\ei{t)) for (jHUl) and = 

■Af±J2k=i^kit)S [Ek{t)]\-&k{t)) for (jn21)- After computing the density operator of these 
cavity- field states, p{t) = |\E' (t)) which refiects all the properties of a quantum sys- 

tem - such as superpositions and decoherence (when fiuctuating parameters are in order) - 
the symmetric ordered characteristic function, defined as in j33|, follows: 

2:5(7,7*,^) = Tr(p(t)e^'^^-^*'^) = {^{t)\e'^^'-^*^\^{t)). (46) 

From the characteristic function ^5(7, 7*, t) we define the Wigner distribution function 

Q 

W{r^, t) = ^^j d'^<ts{l, l\ t) e^*""^^* . (47) 

which will be employed here to represent the quantum properties of the cavity-field state 
conveniently in a three-dimensional Re(?7),Im(?7),iy space. The result of the lengthy and 
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somewhat involved integration over the entire complex plane is presented in Appendix A, 
only for the state (jHUjl . 

Next we analyze the fluctuations of the quadratures of the cavity mode, defined as the 
dimensionless position and momentum operators, Xi = (a + a^) and X2 = —i{a — a^), 
respectively. The dynamic fluctuations for these quadrature operators in the cavity-field 
state, given by 

{AX]) = {X])-{X,)\j = l,2 (48) 
are obtained by computing the variances, 

(Aa^) = (a^) - {af , (49a) 
(A(at)^) = ((at)^)-((at)>\ (49b) 

(Aa^a) = (a^a) - (a^) (a) . (49c) 

We note that with the above definitions for Xi and X2 we get for the coherent state the 
minimum uncertainty value (AX^ = 1. The expected values of the normal ordered opera- 
tors defined in Eqs. ()49a|) . ()49b|) . and ()49c|) may be conveniently evaluated with the help of 
the normal ordered characteristic function 

C:^(7, 7*, t) = Tr (p(t) e^'^' e"^*") = e'^l'/^ ^5(7, 7*, t). (50) 

^From this equation, if we want the normally ordered moments, it is easy to derive the 
expression 

ii-T-n = 5:^^(3^ e:.(7,7*,t)|,.,..o, (51) 



which is suitable for computing the variances in Eqs. (j49aj), ()49b|) . ()49cp . and 



VII. PROTOCOLS FOR THE PREPARATION OF "SCHRODINGER CAT"-LIKE 
STATES AND NUMBER STATES 

A. "Schrodinger cat" -like states 

To prepare a particular superposition state from we follow a three-step protocol. 
1) First, we adjust the amplitude k of the parametric amplification and the atom-field 
interaction time r = t2 — ti in order to obtain a particular angle G = \ifi{t2) — (^2(^2) I defined 
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by the squeezing directions of the states composing the "Schrodinger cat"-hke superposition. 
2) Next, the desired excitation of the prepared state can be achieved by manipulating the 
excitation of the initial coherent state injected into the cavity, the amplitude of the linear 
amplification (that of the parametric amplification has been fixed in the first step), and 
the time interval of the amplification process. 3) Finally, the amplitude of both states 
composing the superposition can be adjusted through the probability amplitudes of the 
atomic superposition state prepared in the first Ramsey zone. 

In Figs. 2(a,b and c) we present some superposition states of the cavity mode generated 
with the above protocol. In all these figures we have considered an atom prepared in Ri in 
the superposition (|0) + |1)) /V2. We have also disregarded the linear amplification process 
while the parametric amplification is switched off at t = t2 when the atom leaves the cavity. 
In the captions of Figs. 2(a,b, and c) we present the fluctuations for the quadrature operators 
and the parameters r(t2) and 9 used for the preparation of the desired states. Figs. 2(a,b) 
indicate the possibility to control the squeezing directions of the quasi-probability distribu- 
tion of the individual states composing the prepared "Schrodinger cat" -like superposition. 
This control will be extremely useful for generating number state as circular squeezed states 
as shown below. 

It is worth noting that a number of exotic reference states have being requested for 



measuring properties of chosen field states. In Ref. the reciprocal-binomial state is 
requested as a reference field for measuring the phase distribution of a chosen field without 
having to obtain sufficient information to reconstruct its complete state. An extension of the 

n 

proposal in Ref. was present for the Q-function measurement where a convenient choice 
of a reference state allows us to measure dispersions of quadrature operators ||5C|| • Therefore, 
we hope that the control of the squeezing directions of the components of superposition states 
achieved through our scheme could also be employed to generate these useful reference states. 
We also mention that the state in Fig. 2(c), considered in the analysis in Section VIII, is 
crucial for achieving long-lived mesoscopic superposition states of the radiation field in cavity 
QED. 

Finally, we recall that the amplification processes could be considered, after the atom- 
field interaction, for controlling the prepared cavity-field state. Both amplification processes 
can furnish excitation to the cavity mode, while the parametric one is able to increase the 
degree of squeezing. 
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B. Number states as circular squeezed states 



/^From the present scheme of atom-driven field interaction it is possible to generate number 
states with higher fidelity than those generated as circular states, i.e., a superposition of M 
coherent states having the same modulus and uniformly distributed around a circle in the 
phase space Isd] . To do this, we have to pass atoms through the cavity, obtaining the state 
defined in Eq. (jHH), where M = 2^. Remembering that Sfc(t) = r^^ ...^^^(t) exp (iyj^^^ ^^^(t)) 
and that the squeezing factor r{t) does not depend on the atomic state, differently from the 
squeezing direction in phase space (pi(t), we get Sjt(t) = r(t) exp {i{pk{t)), where we have 
defined ipk{t) = v9^^^,,,/^(t), with k running from 1 to 2^. With these considerations, we 
obtain from Eq. (|32j) the photon distribution function 

2 



Vnit) = \{n\^^ {t))f 



k=l 



[tanhr(t)]' ^ i[^^(t)-v^{t)]n/2 



' ^' 2"ra!coshr(t) 



k,m= 




\ {\Ut)? + Wm{t)\^) + ^tanhr(t) [{dl{t)f e^^^^'^ + (^^^(t))^-^^-^] | 



2 coshr(t) sinhr(t) j \ a/2 cosh r(t) sinh r(i(:) 



where Hn{x) is the n-th Hermite polynomial evaluated at x. We have assumed that all 
atoms were prepared in the same state (|0) + |1)) / \/2 and detected in their excited states. 

In order to get the superposition of squeezed coherent states centered around the origin 
of the phase space (as required to generate the number state) we have to switch off the 
linear amplification process to obtain, from Eq. ()12b|) . h^{t) = 0, leading to coherent states 
having equal amplitudes 'dk(t) = e ^^''i'---''^n^^^ 0^^ f^('^i) = e"*'^**^*-' a. In addition, assuming 
a is real, we have to adjust Pkit) to (1 — A;)7r/A^ and ipk(t) = {k — l)2n/N, to get a symmetric 
distribution of these states around the center of the phase space. With these adjustments 
(which are achieved through the interaction times between the atoms and the cavity mode 
and also through the parametric amplification parameters K{t),ri{t)) the photon distribution 
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function simplifies to 



Vnit) = Yl 



k,m=l 



[tanhr(t)]" 
n\ coshr(t) 



ap(l-tanh r(t)) 




coshr(t) sinhr(i(:) 



a 



- 2 



(53) 



k,m=l 



In Fig. 3a we present the Wigner distribution function of the state generated from the 
passage of = 2 atoms through a cavity initially prepared in the coherent state |a| = 7.4 
with r = 0.99. These choices of the parameters |a| and r are considered in order to maximize 
the photon distribution function for n = 8, attaining Vn=s{t) = 0.95 which is exactly the 
fidelity |(8| N=2{t))^ of the prepared state with respect to the number state |8). The value 
0.95 is considerably larger than that computed without the parametric amplification process, 
when a circular state is generated with maximum fidelity 0.56 with respect to the number 
state |8). The fidelity 0.56 is computed from an initial coherent state |a| = 2.83. In Fig. 3b 
we plot the Wigner function of the state prepared from the passage of = 3 atoms through 
a cavity initially prepared in the coherent state |a| = 8 with r = 0.67. Here we obtain the 
optimal value Vn=i&{t) = 0.99, to be compared with the fidelity 0.79 computed when the 
amplification process is switched off and a circular state is prepared, from an initial coherent 
state |a| = 3.95. 

Note that with the passage of atoms though the cavity a family of number states is 
obtained: \n = q2^^ with q = 1,2, .... However, we stress that the fidelity of the prepared 
state decreases as the integer q increases. In Table I we present the states \n = q2^) for 
some values of q, in order to compare the maximized fidelities computed from our model 
(jFg = = q2'^\ N=3{t))f) with those derived from the circular states technique (F). We 
do not present the values of |a| and r used to calculated the fidelities. 

We thus conclude that the atom-driven field process is suitable for preparing number 
states with higher fidelity than those generated as circular states. Next we present another 
application of the states generated by the atom-driven field interaction. 

VIII. PREPARING LONG-LIVED MESOSCOPIC SUPERPOSITION STATES 

Evidently, the squeezed superposition in Eq. (jHUj) was ideally prepared. In a real en- 
gineering process the dissipative mechanisms of the cavity and the two-level atom, despite 
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of the fluctuations intrinsic to tlieir interaction, must be taken into account. The complex 
calculations involved in the engineering of quantum states under realistic quantum dissi- 
pation and fluctuation conditions can be computed through the phenomenological-operator 



approach presented in Refs. 



24| . However, in this paper we will not take into account 



the action of the reservoir during the preparation of the squeezed superposition (jSHI)- As 
usual, to estimate the decoherence time, we next consider that an ideally prepared state 
is submitted to the action of a quantum reservoir described by a collection of harmonic 
oscillators whose Hamiltonian is H/j = hukb^b. In addition, we will be interested in the 
action of a vacuum-squeezed reservoir at absolute zero; its initial density operator reads 
Pr = Ylk^k |Ofc) (Ofcl Sl, Sk being the squeezed operator for the k-th bath oscillator mode. 
We are assuming here that, somehow, it is possible to describe all the mechanisms of dissi- 
pation of the cavity in terms of the action of a vacuum-squeezed reservoir. Describing the 
interaction between the reservoir and the system (the cavity mode modeled as = hwa^a) 
as V = ^^,h{\ka)hk + A^afej,), characterized by the strengths , the decoherence time 
deduced from the idempotency defect of the reduced density operator of the cavity field, as 
suggested in [^, is given by 

1^ = - (H)l« + mi)R + mi)s - {^')s,n (54) 

where the Hamiltonian comprehends a sum of three terms H = H5 + B.^ + V. The average 
(H)g ((H)^) is taken with respect to the density matrix of the system (reservoir), given by 
P5 = |\E' (t)) (\I/ (t)\ (pr), where |^ (t)) is given by Eq. From Eq. (jSH), the decoherence 

time of the cavity-field state is given by 

r = , -, , (55) 



{2N + 1) ((at) (a) - (ata)) + 2Re [m ((at)' - ((a^)')) 



N 



where tr is the relaxation time defined by the cavity, = sinh'(r), and M = 
— e*^ sinh(2r)/2, r and (p being the squeeze parameters of the vacuum reservoir j^^. Here the 
mean values are computed from the prepared squeezed superposition (jHDI)- Since the excita- 
tion of the initial coherent state a and the squeeze parameters {r{t2), fe{t2)) have been fixed 
by the engineering protocol, we note that Eq. ()55p depends only on the reservoir squeeze 
parameters {r,ip). Considering the situation where a is real and (a| —a) = exp(— 2a') ~ 
(implying a ^ -\/2), besides the assumption that ^i{t2) = ^ and (^2(^2) = ^ + 2n7i [n integer) 
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(implying = \(pi(t2) — ¥^2(^2)! = ^nn, i.e., the states composing the superposition (I^Uj) are 
squeezed in the same direction), we obtain 

T = tr\1 + cosh(2r) [2a^cosv5sinh(2r) - (l + 2a^) cosh(2r)] 
— sinh(2r) [(1 + 2a^) cos(^ — (^) sinh(2r) — (cos(i^ — 2lp) + cos Lp) cosh(2r) 
+ (cos((^ — 2ip) — cos(^)] I (56) 

where r = r{t2)- The maximization of the decoherence time r in Eq. (j56p with respect to 
the parameters (r,^), leads to the results 

rA=r + + = 0, (57a) 

fs = r - ln(l + W) /4, = TT, (57b) 

which follow when we take v^(t2) = (2m + l)7r and v^(t2) = 2m7r (m integer), respectively. 
When B 7^ 27^7?, the maximum of r turns out to be smaller than that for B = 2n7r, given 
either by the pair (r^, if a) or (re, ^s). Observe that the directions of squeezing of both 
states composing the superposition (jSOI), defined by the angles ^^1(^2) and ip2it2) has to be 
perpendicular to the direction of squeezing of the vacuum reservoir. 

Next, we compute the "distance" in phase space between the centers of the quasi- 
probability distribution of the individual states composing the prepared superposition ()3Up. 
This distance is defined by the quadratures of the cavity field X = (a^ + a)/2 and 
Y = {a — a^)/2i, as 

D = [{{X), - {X),f + ((F), - {Y),f] , (58) 

the subscripts 1,2 referring to the two states composing the superposition. Taking (^1(^2) = 
^^2(^2) = (2m + 1)77 or 2m7r, the distance becomes D = {X)^ — {X)^ = 2Q;exp(r) or 
2aexp(— r), respectively. We will focus on the case ^^1(^2) = ^2(^2) = (2mi + l)n, since it 
results in a large distance D between the two states composing what we actually want to 
be a mesoscopic superposition. The decoherence time and the mean photon number of the 
prepared state, obtained from the values (ja, ^Pa), with exp(— 2a^) ~ 0, read 

r^Tn/a, (59) 
(n) = {^a^a) ^ a^exp(2r) + sinh^r. (60) 



22 



Remarkably, with the approximation exp(— 2a^) ^ 0, the decoherence time for the prepared 
cavity-field state when v^i(t2) = V^2(^2) = (2m2 + l)7r — under the action of a vacuum 
reservoir squeezed in the direction if^ = — turns out to be practically independent of 
the parameter r and thus of its intensity (n) and distance D. Therefore, the decoherence 
time (|59p becomes practically independent of the quantities which define the macroscopic 
character of the cavity-field state. From the result in Eqs. (jSHl) and (plj) we conclude 
that it is convenient to start from a coherent state a as small as possible (within the limit 
exp(— 2«^) ^ 0) and to adjust the macroscopic coupling parameter l^^l in order to obtain 
a large squeeze factor and, consequently, a large intensity of the prepared state and a large 
distance D, since we are actually interested in mesoscopic superpositions. We stress that 
even considering the weak coupling regime (|^^| < 1) we obtain, from Eqs. (j36a|) and (j42j) . 
large squeeze parameters: considering \Pi\ = 0.1, a = ^/2, and the experimental running 
time about 2 x 10~^s, we get a superposition state where r ^ 2 and (n) ^ 10^ photons. 

The mechanism behind this result is the degree of entanglement between the prepared 
state and the modes of the reservoir, which depends on the relative direction of their squeez- 
ing, defined by the angles ^^1(^2) = V^2(^2) and (fA- A result supporting this argument is 
presented in j52| where it is shown that the injection of two modes, squeezed in perpendic- 
ular directions, in a 50/50 beam splitter does not generate an entangled state. A careful 
analysis of the dependence of the relative direction of squeezing on the degree of entangle- 
ment between a prepared state and its multimode reservoir will be presented in [5^ . Despite 
the fact that the mechanism behind the long-lived mesoscopic superpositions is mainly the 
perpendicular squeezing directions between the prepared state and the reservoir modes, the 
magnitude of the parameter r plays a crucial role in the present scheme for producing the 
mesoscopic superposition by increasing both its intensity (n) and distance D in phase space. 

The values presented above for r, (n), and D are to be compared with those when 
considering a non-squeezed (NS) cavity-field state {{n)j^g = a^, D]\rs = 2a ) under the 
infiuence of i) a squeezed reservoir, resulting in the decoherence time Tj ~ T/j/a, and ii) a 
non-squeezed reservoir, such that th ^ Tji/2a'^. Note that in both cases i) and ii) we obtain 
the ratios {n) / {n)j^g ^ exp(2r) and D/Dns ~ exp(r). Therefore, despite the exponential 
increase in the ratios of both excitation and distance, we still get r ^ Ti when comparing 
our results with previous schemes in the literature, where a squeezed reservoir is assumed 
for the enhancement of the decoherence time for non-squeezed cavity-field states and 
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reservoir, we obtain a still better result r ^ aru. 

IX. DISCUSSION AND CONCLUSION 

We have presented a scheme for the preparation and control of a cavity-field state through 
atom-driven field interaction. The Lewis and Riesenfield time- dependent invariants |48| 
were employed to obtain the eigenstates of the cavity mode dispersively interacting with 
a two-level atom and simultaneously under linear and parametric amplification processes. 
Protocols for preparing particular superposition states and the number state were presented. 
While relying on the manipulation of the initial states of the cavity mode and the two-level 
atom, considered in previous schemes [s^, our protocol also employs the time-dependent 
parameters involved in the amplification sources to achieve particular superposition states 
and number states. We plotted some interesting "Schrodinger cat" -like states and number 
states generated as circular squeezed states. We demonstrated that the number states 
generated as circular squeezed states exhibit higher fidelities than those generated as circular 
states. 

We have shown how to prepare truly mesoscopic "Schrodinger cat" -like states of the cavity 
field, actually squeezed superposition states, through their coupling to likewise squeezed 
reservoirs. When assuming that the squeezing direction of the cavity field is perpendicular to 
that of the reservoir modes, we found that the decoherence time of the prepared superposition 
state depends only on the initial coherent state of the cavity field from which the squeezed 
superposition is generated. Therefore, the decoherence time is independent of the average 
photon number and the distance in phase space between the centers of the quasi-probability 
distribution of the individual states composing the squeezed superposition. This result 
follows from the degree of entanglement between the prepared state and the modes of the 
reservoir, which depends on the relative direction of their squeezing. When the squeezing 
direction of the prepared superposition and that of the reservoir modes is perpendicular, the 
noise injected from the reservoir into the prepared cavity mode is minimized. A detailed 
analysis of the dependence of the relative direction of squeezing on the degree of entanglement 
between a prepared state and its multimode reservoir will be presented in 0]. 

The experimental implementation of the proposed scheme relies on the possibility of en- 
gineering a squeezed reservoir as well as of parametrically driving cavity-field radiation. We 
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stress that a scheme to reahze physically a squeezed bath for cavity modes, via quantum- 
nondemolition-mediated feedback, has already been presented in Ref. j56|. However, the 
feedback process in does not eliminate the standard nonsqueezed bath and, as we have 
pointed out, our scheme requires an optimal squeezed-vacuum reservoir. The subject of 
quantum-reservoir engineering has attracted some attention, specially in the domain of 
trapped ions 2^; more specifically, a scheme has been presented for the engineering of 
squeezed-bath-type interactions to protect a two-level system against decoherence j^]- 

We emphasize that a proposal to implement the parametric amplification of an arbitrary 
radiation- field state previously prepared in a high-Q cavity is presented in Ref. 0]. As 
mentioned above, in this proposal the nonlinear process is accomplished through the dis- 
persive interactions of a single three-level atom simultaneously with a classical driving field 
and a previously prepared cavity mode whose state we wish to squeeze. Moreover, regard- 
ing parametric amplification of cavity fields, a technique was recently suggested, based on 
pulsed excitation of semiconductor layers (on the cavity walls) by laser radiation [^. It is 
worth mention that all the treatment developed above in the context of cavity quantum elec- 
trodynamics, for delaying the decoherence process of a squeezed superposition by coupling 
it to a vacuum-squeezed reservoir, can also be implemented in ion traps. We finally men- 
tion that the proposal presented here should motivate future theoretical and experimental 
investigations. 

A. Appendix A 

In this appendix we present the Wigner function computed from Eq. ()47|) and the relation 

Csij, 7*, t) = el^lV^ C:^(7, 7*, t) = el^l'/2 (t) | e"^*'^ e^'^' \^ (t)), (61) 
derived from the antinormal ordered characteristic function 

(2:^(7,7*,t) = Tr(p(t)e-^*»e^'^') = J ^ v]/ (t)) |%-7Wt (62) 

First we have to compute, from Eq. (jSOI), the final cavity- field state I^E' (t)) which, after a 
lengthy calculation, becomes 

\^ (t)) = \^ (t)) = ^^+ [ciri(t)5(£i(t)) \9^it)) + c2T2it)s{e2it)) le^m (63) 
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where the TD function Ti{t) is defined in terms of that given by Eq. (fT^ . as 

r,(t) = T,(t)T,(t2)T,(ti), (64) 

and Oiit) is defined by Eq. (fT3j). From Eqs. (jU^ and we obtain, with i,j = 1,2, the 
result 



X exp 



{ai + 7)(a* + 7*) + bi{a* - + b*{ai + -ff 



1 - ^bib* 



(65) 



where the TD function Kj^ reads 



K,^ = \Af+f c,c*T,{t)Tm sech(r(t)) 



X exp 



|2 



and 



-1 (l^.(t)l^ + \eM") + ^tanhr(t) (e^-^W (^*(t))^ + e'^^^W {9,{t)f) 



di = Oi(t) sechr(t), 

hi = -^tanhr(t) e'^^W 



(66) 

(67) 
(68) 



Note that for the weak couphng regime, Ti(t) = r2{t) = r{t). Finally, from Eq. (j47|) and the 
characteristic function in Eq. ()65|) we obtain the Wigner function 



exp 



where the TD functions Aij, Bij, Cij, Dij,.Eij, and Fij satisfy 



(69) 



A 



1 - Abib: 



exp 



ai{t)a*{t) + ai{t)b*{t) + {a*{t)f hit) 
l-Ah{t)b*At) 



B 



Cij 



1 



2 l-46,(t)6*(t)' 
6.(t) 



D 



l-Ab,{t)b*{t)' 

m 



l-4bi{t)b*{t)' 

2ai{t)b*{t) + a*{t) 



-77 + 



l-46,(t)6*(t) 



r/ - 



2a*{tMt) + a,{t) 
1 - Abiitmt) 



(70) 
(71) 
(72) 
(73) 
(74) 
(75) 
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Figure Captions 

Fig.l. Sketch of the experimental setup for atom-driven field interaction. 

Fig. 2a. Wigner function obtained when Q — a = 0, k — x/20, and x't = 2.06. 
The variances for the quadrature operators read (AXi) = (AX2) = 2.63 and the squeezing 
factor attained is r{t2) = 1.45. The mean photon number is 2.96. 

Fig. 2b. Wigner function obtained when Q — n, a — k — x/20, and xr = 3. 99. The 
variances for the quadrature operators read (AXi) = 5.58 and (AX2) = 5.93 and the 
squeezing factor is r{t2) — 1.55. The mean photon number is 26.38. 

Fig. 2c. Wigner function obtained when Q — 0^ a — 5, k, — x/20, and x^ — 1-68. 
The variances for the quadrature operators read (AXi) = 0.31 and (AX|) = 32.0 and the 
squeezing factor is r{t2) = 1.5. The mean photon number is 256.17. 

Fig3a. Wigner function of the state generated by passing N — 2 atoms through a cavity 
initially prepared in the coherent state \a\ — 7.4 with r = 0.99, leading to the maximized 
photon distribution function Vn=?,{t) — 0.95. 

Fig3b. Wigner function of the state generated by passing N — 3 atoms through a cavity 
initially prepared in the coherent state |q;| = 8 with r = 0.67, leading to the maximized 
photon distribution function Pn=i6(^) = 0.99. 
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Tables 

Table I. The fidelity of the states |^'Ar=3(^)) generated from our model {Tq) and those 
derived from the standard circular states technique (F), for different values of the desired 
number states \q2'^). 

q T F 
1 0.99 0.98 
3 0.98 0.65 
5 0.96 0.50 
10 0.90 0.35 
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